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ABSTRACT 

An anomalous radio continuum component at cm-wavelengths has been observed in 
^3 various sources, including dark clouds. This continuum component represents a new 

property of the ISM. In this work we focus on one particular dark cloud, the bright 
reflection nebula M 78. The main goal of this work is to investigate the cm-wave 
continuum emission in a prominent molecular cloud, nearby and with complementary 
observational data. We acquired Cosmic Background Imager (CBI) visibility data of 
M 78 at 31 GHz with an angular resolution of ~ 5.8', and CBI2 data at an angular 
resolution of ~ 4.2'. A morphological analysis was undertaken to search for possible 
correlations with templates that trace different emission mechanisms. Using data from 
WMAP and the Rhodes/HartRAO 2326 MHz survey we constructed the spectral en- 
ergy distribution (SED) of M 78 in a 45' circular aperture. We used results from the 
literature to constrain the physical conditions and the stellar content. The 5 GHz - 
I 31 GHz spectral index in flux density (a = 1.89 ± 0.15) is significantly different from 

• • optically thin-free-free values. We also find closer morphological agreement with IR 

m ^ dust tracers than with free-free sources. Dust-correlated cm-wave emission that is not 

k> due to free- free is significant at small scales (CBI resolutions). However, a free- free 

j_j background dominates at cm- wavelengths on large scales (~ 1 deg). We correct for 

C$ this uniform background by differencing against a set of reference fields. The differ- 

enced SED of M 78 shows excess emission at 10-70 GHz over free-free and a modified 
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blackbody, at 3.4 a. The excess is matched by the spinning dust model from Draine 



& Lazarian ( 1998 ) 



Key words: radio continuum: ISM, radiation mechanisms: general, infrared: ISM, 
ISM: dust, ISM: clouds 



1 INTRODUCTION 

Since 1996, experiments designed to measure the CMB 
anisotropy have reported an anomalous diffuse foreground 
in the range of 10-60 GHz. This diffuse emission is corre- 
lated with thermal emission from dust grains at 100 /im. The 
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spectral index (considering S v oc v a ) of the radio-IR corre- 
lated signal is a ra dio/iR ~ in the range 15-30 GHz, as for 
optically thin free- free ( Kogut et al.|1996 l. But a radio/nt ~ 
—0.85 between 20-40 GHz, for high-latitude cirrus (Davies et 
al.|2006 \ . Additionally the absence of Ha emission concomi- 
tant to radio free-free emission, would require an electron 
temperature T e ^ 10 6 to quench Hi recombination lines 
(Leitch et al 19971. Another emission mechanism was pre- 
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sented by Draine & Lazarian (19981, who calculated that 



spinning interstellar dust grains produce rotational dipole 
emission in the range from 10 to 100 GHz, at levels com- 
parable to the excess diffuse foreground observed over the 
free-free component. 

Observations of specific targets may shed light on 
the anomalous foreground, whose existence in the dif- 
fuse ISM has been inferred statistically. Anomalous cm- 
wavelength radiation lias been found in about a dozen clouds 
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The search by |Finkbeiner et al ( 2002 1 for anomalous 



cm-wave emission in 9 dark clouds, resulted in one detec- 
tion: LDN 16223 The radio spec trum of LDN 1622 matches 



that of spinning dust emission (Finkbeincr 2004). Cosmic 
Background Imager (CBI, |Padin et al.||2002[ ) observations 
of LDN 1622 linked, on morphological grounds, the cm- 
wave emitters to the smallest interstellar dust grains (i.e. 
the VSGs) - as required for spinning dust (' Casassus et al. 



20061. Interestingly, LDN 1622 is notable in being exposed 
to UV radiation from the adjacent hot stars in the Ori OB lb 
association. 

The prototypical dark cloud LDN 1688, close to the 
p Oph star, was recently found to be bright at cm- 
wavelengths (ICasassus et al. 20081. The p Oph cloud is 



undergoing intermediate-mass star formation. UV radiation 
from its hottest young stars heats and dissociates exposed 
layers, but does not ionize hydrogen. Only faint radiation 
from the Rayleigh- Jeans tail of ~ 10-100 K dust is expected 
at wavelengths longwards of ~3 mm. Spinning dust com- 
fortably explains the radio spectrum of p Oph W, as in 
LDN 1622. However, spinning dust encounters difficulties 
in explaining the cm-wave morphology of the p Oph cloud. 
Diffuse mid-IR (10-20 /im) emission from p Oph is inter- 
preted as stochastic heating of VSGs (Bernard et alll 1993 1 . 



Dust emissivities are proportional to both the dust grain 
density and the local UV energy density. Spinning dust was 
introduced to explain the 100 fim-1 cm correlation in the dif- 
fuse ISM, but does this correlation extend to smaller scales 



in denser environments? Casassus et al. (20081 report that 



the mid-IR intensity peaks in p Oph, from the circumstel- 
lar nebulae about Oph S 1 and Oph SR 3, have no 31 GHz 
counterparts. They find that the predicted spinning dust in- 
tensities towards S 1 are in excess over the observed values 
by a factor of > 40. On the basis of Spitzer IRS infrared spec- 
troscopy, this discrepancy could marginally be explained by 
VSG depletion near S 1 and SR 3. 

As part of an effort to understand the emission mecha- 
nisms in environments giving rise to the 31 GHz continuum, 
we have studied a dark cloud with extensive complemen- 
tary data: the reflection nebula M 78, with good constraints 
on its physical environment. Several observations are avail- 
able in the literature that constrain the physical conditions 
in M 78, making it a promising testbed of candidate emis- 
sion mechanisms. As a reflection nebula M 78 has relatively 
low levels of free- free emission but also contains stars bright 
enough to excite circumstellar dust grains. Reflection neb- 
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Figure 1. M 78 area taken from DSS blue band Survey. The 
largest and brightest nebula is NGC 2068/M 78, the arc-shaped 
nebula to the west is NGC 2067. Northeast from NGC 2068 is 
NGC 2071 and the faintest nebula, southwest from NGC 2068 
is NGC 2064. The crosses indicate the positions of the bright- 
est stars in the field. The dashed lines point to these stars and 
the nebulae are pointed to with solid lines. The coordinates are 
in right ascension and declination offsets from (05 ,1 46 m 46 s .7 + 
00°00'50" J2000), in degrees of arc. 



ulae are regions of star formation, and they contain large 
amounts of dust. 

In this report we refer to M 78 as a group of reflec- 
tion nebulae, all of them part of LDN 1630. This is a large 
dark cloud in the constellation of Orion, north from the belt, 
which also contains NGC 2023. LDN 1630 corresponds to the 
Orion B molecular cloud, which has recently been studied 
in molecular lines by |Buckle et al.| ( |2009[ ). The area that 
hereafter will be called M 78 contains the following nebulae: 
NGC 2068 (sometimes also referred to as M 78), NGC 2071, 
NGC 2064 and NGC 2067 (see Figure [lj. 

Section [2] refers to the physical conditions and different 
phases of gas and dust found in M 78. In Sec. [3] we will 
describe the observations used for this work. Sec. [3] also in- 
cludes a discussion about the CBI flux loss. In Sec. [4] we 
present our results, derived from a comparison between the 
CBI data and various tracers of dust and free-free emis- 
sion. In Sec. [5] we summarize our conclusions, quantifying 
the 31 GHz radiation from M 78. 



2 ENVIRONMENT AND PHYSICAL 
CONDITIONS 

L1630 is embedded in the Orion Complex, at a distance of 



about 400 pc ( Anthony- Twarog 19821. It is part of Orion 



Molecular Cloud II, among other star forming regions and 
recently formed clusters. We also find Barnard's Loop in 
this region; it passes at ~ 40' from the northern end of 
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NGC 2071. As a reflection nebula M 78 has low levels of free- 
free emission but is exposed to stars bright enough to excite 
IR emission from dust. Reflection nebulae are regions of star 
formation, so they contain young stars and dense molecular 
cores. As we can see in Fig. [2j although the large scale radio 
emission is dominated by a uniform free- free background, the 
individual features of M78 become more pronounced with 
increasing frequency and angular resolution. 

The average H- nucleus density value in LI 630, derived 
from CO low resolution maps, is n ~ 10 3 cm -3 , but to- 
wards NGC 2068 and NGC 2071, higher densities are found 



(Strom et al. 19751. Lada et al. (19971 obtained, for these 



clumps, densities of log(n/cm 3 ) = 5.3 for NGC 2068 and 
log(n/cm~ 3 ) = 5.9 for NGC 2071. These values were de- 
rived from an excitation analysis of CS line fluxes for three 
transitions J = 2 — > 1, J = 3 — > 2 and J = 5 — > 4. Buckle et 
al. | |2009| ) have recently reported on a large scale CO(3-2) 
isotopologue survey of Orion B. including the M 78 area, 
and focusing on the star formation activity. Depending on 
the isotopologue, the molecular mass of NGC 2071 is 0.4 to 
3.7 10 3 Mq. NGC 2071 dominates NGC 2068 b y mass. 

M 78 is still undergoing star formation (Flaherty & 



Muzerolle 20081. M 78 contains two embedded clusters, in 
both NGC 2071 and NGC 2068. These clusters are consid- 
ered to have an age of about 2± 1 Myr, and for a significant 
fraction of stars in the cluster there is evidence of accretion 
disks. 

The illuminating stars for NGC 2068 and NGC 2071 



have been described by Strom et al. ( 1975 1. NGC 2068 hosts 



the star HD38563N (B2 m) and embedded in NGC 2071 is 
HDE290861 (B2-B3). Other early type stars in NGC 2068 
are HD28563S (B3-B5) and HD38563C (AOll). These stars 
appear to be very young, with ages around 10 5 years. For 



NGC 2064 the illuminating star is SSCV 111 (Strom et al 
1975[ ), and its spectral type is B3v (Chini et al.|1984[ ). 

The photo-dissociation-region (PDR) surrounding 



HD38563N was studied by |Owl et al.| ( |2002[ ), finding a dust 
temperature of 41 K, a gas density of 5000 cm -3 and a gas 
temperature of 250 K. They also give Go = 2800 for the UV 
radiation incident on NGC 2068, where G a = 1 corresponds 
to the ave rage interstellar ra diation field of 1.6 10 -3 erg 
cm" 2 s" 1 (Mathis et al.|l983 l. 

M 78 has an associated C n region. According to 
Pankonin & Walmsley ( 1978 1 the C II region extends ~ 20' 



along a line between NGC 2071 and NGC 2068. The elec- 
tron temperature ranges fro m 20 K to 50 K and the elec- 
tron density is 0.2-1.0 cm -3 ( Brown et al.|[l975 Pankonir 
|fc Walmsley|[l978l ). Both densities and temperatures were 
derived from the ratio of the power in two carbon recom- 



bination lines. According to Silverglate & Peter ( 1983 ) this 



C II region is unusual because it does not present a sulfur 
line in its spectrum. This sulfur line is found towards all 
other known C II regions. Another datum is that the carbon 
radio recombination line is the narrowest yet observed for 
reflection nebulae and has fainter brightness temperatures 
compared to other Cn regions. 



Besides the Cn regions, a compact Hn region (~ 1') 


was reported by[Matsakis et al. 


(|1976|) around the brightest 


star in NGC 2068, HD38563N. 


Matsakis et al. ( 1976 1 mea- 



sured 0.1 Jy at 2.4 GHz. By assuming a temperature of ~ 
7000 K, they found an electron density of 40 cm -3 . This 
region is responsible for most of the free-free emission of the 



nebula on 1-10' scales, although on 1 deg scale it is weaker 
than the background diffuse H II region. 



M 78 also hosts H 2 masers ( |Campbell||1978[ ). These 
masers are considered tracers of proto-stellar disks, given 
that they coincide with infrared sources with associated IR 
excesses (i.e. circumstellar disks) and collimated outflows 
(e.g. as in class-II protostellar objects). 



3 OBSERVATIONS 

3.1 Cosmic Background Imager data 



The CBI is described in Padin et al. (2002). The CBI is a 



planar interferometer array with 13 antennas, each 0.9 m 
in diameter, mounted on a 6 m tracking platform. The 
CBI receivers operate in 10 frequency channels covering 26- 
36 GHz. It is located in Llano de Chajnantor, Atacama, 
Chile. Cancellation of ground and Moon contamination was 
obtained by differencing with a reference field at the same 
declination but offset in hour angle by the duration of the 
on-source integration. We used an on-source integration time 
of 8 min, with a trailing reference field. 

We used the CBI to observe M 78 on 20 December 
2004. The primary beam is 45'. 2 at FWHM of the best fit- 
ting Gaussian at the reference frequency of 31 GHz, and the 
synthesized beam is 5'. 94 x 5'. 70 with uniform weights. We 
also observed M 78 with CBI2, an upgrade where the 0.9 m 
CBI dishes were changed for 1.4 m dishes. In the case of 
CBI2, since the primary beam FWHM is 28'. 2 (this beam 
diverges from a Gaussian outside ~FWHM/2), we mosaiced 
two pointings on NGC 2071 and NGC 2068. The uv cover- 
age was different in the two fields, thus they have different 
synthesized beams (using uniform weights); for NGC 2068 
we have 4'. 89 x 3'. 97, while for NGC 2071 the values are 
4'. 37 x 4'. 00. The visibility coverage in the uv plane is dis- 
played in Fig. [3] 

The data were reduced and edited using a special- 
purpose package (CBICAL, developed by T.J. Pearson). 
Flux calibration was performed using Tau A, whose fluxes 
are in turn calibrated against Jupiter (with a temperature of 
146±0.75 K,|Hillet al. 120091). The flux calibrator is also used 



as the reference for an initial phase calibration (restricting 
to baselines longer than 2 m to avoid ground spill over). 

The electronic phase stability of the CBI system is bet- 
ter than 10 deg over the timescales between primary phase 
calibrations. However the pointing accuracy of the mount is 
approximately 0.5 arcmin root-mean-square (rms). Due to 
the co-mounted nature of the CBI antennas this appears as 
correlated phase errors in the visibilities and translates to 
an equivalent pointing uncertainty in the map. The pointing 
error can be measured, and removed, by interleaved observa- 
tions of a nearby phase calibrator. For the CBI observations 
on 2004 Dec 20, the calibrator J0541-051 was observed, and 
the resulting pointing error correction was 62 arcsec. For the 
CBI2 observations no sufficiently bright calibrator was ob- 
served. We note that these pointing errors and corrections 
are in any case a small fraction of the synthesised beam size 
(6 arcmin for CBI and 4.5 arcmin for CBI2) and do not sig- 
nificantly affect the positional matching of features in the 
CBI maps with data at other wavelengths described below 
(see |Casassus et al.|2 006 for tests on the pointing accuracy 
of a similar dataset). 
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Figure 2. Large scale structure around M 78 in various frequencies. The solid line circle is the aperture where we integrated to extract 
flux densities. The crosses indicate the center of the reference fields used for differentiation. The bright source at ~ 2.5 degrees southwest 
from M 78, common to each image, is NGC 2023. (a) At 2326 MHz, from Rhodes/HartRAO survey, we see clearly the arc to the northeast 
which corresponds to part of the Barnard's Loop, but there is no structure at the position of M 78. (b) At 93.4 GHz, corresponding 
to WMAP W band, we see a faint feature corresponding to M 78. (c) Image at 100 fim from IRIS, we see that the structure of M 78 
corresponds to that observed in (b). x— and y— axis show right-ascension and declination offsets from (05 ,l 46 m 50 s .4 + 00°09'31" J2000), 
in degrees of arc in degrees of arc. 
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Figure 3. The uv coverage for (a) CBI when taking M 78, (b) for the observation of NGC 2068 with CBI2 and (c) for CBI2 in NGC 2071. 



We fitted the CBI visibility data with 3 elliptical 
Gaussians, corresponding to NGC 2068, NGC 2071 and 
NGC 2064. The restored image is shown in Fig. [4^. In the 
case of CBI2 the components are coincident with the four 
nebulae NGC 2068, NGC 2071, NGC 2067 and NGC 2064. 
The Gaussian model convolved with the synthesized beam 
plus residuals of the CBI and CBI2 data are shown in Fig. [4^. 
and Fig. [5^,, respectively. Hereafter we refer to these images 
as "model-fitted" images. The CBI2 image shows clearly the 
improved resolution of CBI2, and the increased flux loss of 
CBI2 due to reduced coverage at short spacings. 

For CBI2 we also implemented the MEM algorithm (as 
in |Casassus et al.|2006[ ). The MEM reconstruction provides 
a model sky consistent with the data, prior to convolution 



with the synthetic beam. We show a comparison of the two 
models in Fig. [5] 



The expected thermal noise for the CBI image is 12.78 
mjy beam -1 , for CBI2 we have two values, one for M78 
of 9.94 mjy beam" 1 and another for NGC2071, of 11.91 
mjy beam -1 . The achieved noise can be measured by taking 
the rms dispersion of a region within the half-power contour 
of the mosaiced primary beams, before correcting for it. The 
noise in Fig. [4^,, for CBI1, is 15 mjy beam -1 , and the noise 
in Fig. [H^,, for CBI2, is 20-30 mjy beam -1 (the combina- 
tion of two fields results in a varying noise map, we take 30 
mjy beam -1 as a conservative value). 
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Figure 4. Comparison of CBI contours and various templates. The units are Jy beam -1 . The dashed circle corresponds to the CBI 
primary beam, the small circle outlined is the CBI synthesized beam, (a): restored CBI image. Contours are taken at 95, 80, 65, 52, 43 
and 35% of the peak emission. White crosses indicate the brightest stars: from north to south, HDE290861, HD38563N, HD38563S and 
HD38563C. (b), (c), (d) and (e): restorations of CBI-simulated visibilities on the IRIS images at 12, 25, 60 and 100 fim, with contours 
from CBI. (f): PMN image at 4.85 GHz with CBI contours. Coordinates follow from Fig. [5] 



3.2 Auxiliary Data 

The images that are used to make a morphological compar- 
ison were taken from IRIS ( jMiville-Deschenes fc Lagache 
20051, a re-processing of the IRAS survey (Wheelock et al 



1991 1 at 12, 25, 60 and 100 fim. We use these infrared data 



to test for correlation between the 31 GHz emission and 
the emission from interstellar dust. The mid-IR emission at 
12 |im is thought to be produced by stochastic heating of 
very small grains (VSGs; e.g. |Draine fc Li| ( |200l] )). The heat 
capacity of these grains is small enough so that the absorp- 
tion of a single UV photon will increase its temperature to 
above 100 K, but their small cross section makes it unlikely 
that they encounter many simultaneously. Thus VSGs are 
not in equilibrium with the radiation field. At longer IR 
wavelengths, (such as in the 60 and 100 fim IRIS bands), 
the emission is due to classical grains at T ~ 20 — 50 K. 



The resolution of the IRIS data is not fine enough to 
allow a comparison with CBI^] To trace small-scale struc- 
ture, that could be compared with the improved resolution 
of CBI2, we used archival data acquired with the IRAC in- 
strument aboard the Spitzer Space Telescope. This image 
is taken at 8 fim and traces VSGs and PAHs, like IRIS at 
12 fim, but with a finer resolution of 4" that allows proper 
point-source subtraction and comparison with CBI2 data. 
The IRAC image is shown in Fig. [6^,. 

Point sources from IRAC were subtracted using a me- 
dian filter. For the remaining sources, they were subtracted 
individually by fitting a Gaussian to the source and then 



2 as the IRIS resolution, of 4 arcmin, approaches that of CBI2, 
the simulation of CBI2 observations on the IRIS images will yield 
coarser resolutions than real CBI2 maps 
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Figure 5. Comparison of model-fit and MEM reconstructions of CBI2 data for M 78. Contours were taken at 95, 80, 65, 52, 43 and 
35% of the peak emission for each image. The dashed line corresponds to a noise level of twice the minimum noise, (a) Model-fitted 
reconstruction with contours. The maximum in the image is 0.32 Jy beam -1 , (b) MEM model with contours. The maximum is of 
0.23 MJy sr -1 . (c) Restored MEM image, the model was convolved with the clean beam and residuals were added, the peak intensity 
in this image corresponds to 0.32 Jy beam -1 . All images are corrected for the primary-beam (including the MEM model). Coordinates 
follow from Fig. [2] 



subtracting it. CBI observations were simulated on both the 
IRIS and IRAC images using the MOCKCBI package (Pear- 
son 2000, private communication), in order to have equiva- 
lent visibility sampling when comparing with CBI and CBI2 
images. This will remove any uniform background from IRIS 
and IRAC. When simulating IRAC data we embedded the 
image in a background taken from IRIS 12 fim, to extend 
the field of the IRAC 8 (im mosaic to at least 3 times the 
size of the primary beam. When embedding IRAC 8 fim 
in IRIS 12 fjja we fitted a plane to the difference of emis- 
sion from the border of the images in order to remove most 
of the discontinuities and to remove the contribution from 
zodiacal light in the 8 (im image (IRIS images have only 
residual contribution from zodiacal light). The simulation is 
also useful to acquire visibility-plane data, which we used 
to make linear correlations to obtain a quantitative compar- 
ison between dust and 31 GHz emission. The simulations 
obtained will hereafter be referred to as CBI-simulated or 
CBI2-simulated images. 

We also carried out a morphological comparison of the 
CBI2 data with PMN (at 4.85 GHz) flCondon et al.||1991[ 



survey at 2.326 GHz with a resolution of 20' (Jonas et al 



with a resolution o f 4.1', see Sec. |4.6|), NVSS (1.4 GHz) 
( Condon et al.|1998| with a resolution of 45") and SHASSA 



Gaustad et al. 



2001) 



(at 6563A, and with a resolution of 48 
These templates trace free-free emission, which is present, 
at some level, at 30 GHz. It is necessary to quantify the 
intensity and spatial distribution of the free-free component 
in order to test for the existence of an excess at 30 GHz. 
These comparisons are discussed in Sec. |4.6| 

For measurements of flux density we have used images 
taken from WMAP which give information about the con- 
tinuum radiation on larger scales. The WMAP channels are 
at 22.8, 33.0, 40.7, 60.7 and 93.5 GHz, with resolut ions of 
0.88, 0.66, 0.51, 0.35 and 0.22 degree respectively ( |Limon| 
et al.||2009| . We also used data from the Rhodes/HartRAO 



19981 



3.3 Flux Loss Correction and Spectral Energy 
Distribution 

The CBI, being an interferometer with no total power mea- 
surements, is insensitive to an uniform background. The lack 
of total power also removes part of the extended emission, 
thus high-pass filtering the sky intensity field. 

In order to estimate the level of flux loss, due to missing 
short spacings, we used the IRIS images as templates, shown 
in Fig. [4] to correlate with the CBI signal. We calculated 
the integrated fluxes in a 45' aperture from both raw and 
CBI-simulated data at each of the IRIS wavelengths, and 
compared them. This showed an average flux loss of 34 ±8% 
among the IRIS bands. This means that most of the flux in 
the IRIS images is within the compact objects corresponding 
to NGC 2071, NGC 2068 and NGC 2064. 

Table [1] summarises the SED of M 78, for flux densities 
extracted in a circular aperture 45' in diameter, correspond- 
ing to the CBI primary beam diameter, and centered on RA: 
05' 1 46 m 46 s .7 and DEC: +00°00'50". Our choice of aperture 
contains the bulk of the M 78 region, and matches the signal 
seen by CBI. Using the same aperture diameter we also cal- 
culate the flux density in 8 reference fields which are marked 
as crosses in Fig[2] Subtracting the flux density in each refer- 
ence field from the flux density centered on M 78, we remove 
the uniform background, and estimate the error in this back- 
ground level with the root-mean-square (rms) dispersion of 
the 8 different background-subtracted flux densities. 

The error bars in Table [I] give 1 a, with a — 



2 

noise 



+ cr| iff , where: 



• Odiff is the rms dispersion of the results from differenti- 
ations against the different reference fields, 
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Table 1. Flux densities extracted in a circular aperture 45' in 
diameter, and centred on M 78. 



v a 


F i 




2.326 


7.08 ±0.21 


0.89 ±0.48 


22.8 


4.40 ±0.31 


2.34 ± 0.66 


33.0 


4.34 ± 0.47 


2.79 ±0.81 


40.7 


4.37 ±0.59 


3.18 ±0.98 


60.7 


5.07 ±0.69 


4.68 ± 1.44 


93.5 


20.20 ± 1.47 


19.55 ± 3.60 


31 


3.96 ± 0.30 d 


2.73±0.28 d 


3000 


32418 ± 3220 


29700 ± 3310 



"Frequency of observation in GHz. 'Flux density in Jy. c Flux 
density in Jy after differencing with the reference fields. d This 
value is obtained from IRIS 100 /mi-31 GHz correlated emission, 
using the slope found from the visibility correlation. e Origin of 
the data, for each frequency (in GHz): 2. 326, |Jonas et al7| ( |1998] ) ; 
22.8, 33.0, 40.7, 60.7, 93.5,|Limon et al.|<|2009| l; 31, this work; 
3000, |Miville-Deschenes fc Lagache| | |2005| i 



• ^noisc = Af2/ rms y'iVbeam Alport , where An is the solid 
angle of one pixel, 7 rms is the nominal thermal noise in spe- 
cific intensity, when available, or the root-mean-square dis- 
persion of a ~0.5 deg square box devoid of signal, iVbeam is 
the number of pixels in a synthesized beam for each image, 
A^apert is the number of pixels in the aperture, in this case a 
diameter of 45'. 

An additional term <j ca i could have been added in quadra- 
ture to the 1 a errors on a flux density F, <7 ca i = 0.1 F, as a 
representation of systematic errors in calibration. However 
tTdiff dominates the error budget, and o" ca i is negligible when 
added in quadrature. 



4 DISCUSSION 

To characterize the emission from M 78 that may not 
be explained as free-free and/or vibrational dust, i.e. the 
Rayleigh-Jeans tail of sub-mm dust emission, we will per- 
form a morphological and spectral analysis. First, we will 
summarize a comparison in the sky plane with various tem- 
plates that trace different phases of dust and gas, to find 
qualitatively which corresponds better with the emission 
found at 31 GHz. Secondly, we will correlate the CBI data 
with the CBI-simulated templates, both in the visibility and 
sky planes. The visibility cross-correlations give informa- 
tion on the average properties over the CBI primary beam. 
Thirdly, another diagnostic is fitting a SED. This will give 
us the spectrum of M 78, which allows us to estimate free- 
free and vibrational dust contributions at 31 GHz, and infer 
the level of any excess. 



4.1 Qualitative CBI - IR comparisons 

In Fig. |4j3-e we compare the CBI and IRIS images. There is 
a qualitative correlation between CBI and IRIS, with some 
variations. At 25 (im NGC 2071 is more prominent than 
NGC 2068, which contrasts with IRIS 12, 60 and 100 /mi. 
For the other IRIS images it is difficult to tell by visual 
inspection which correlates best with CBI. 

In Fig. [6] we compare the CBI2 MEM model in contours 



with the IRAC 8 /im image in grey scale. This shows us that 
the spatial correlation obtained between IRIS 12 /jm and 
CBI still holds on smaller scales. We also see a radio peak 
with no IR counterpart towards NGC 2067. In IRAC this 
nebula appears as diffuse emission, while in CBI2 this peak 
is more intense than the one observed towards NGC 2064. 
In Fig. [6jj we show the reconstructed CBI2-simulated image 
from IRAC obtained by model-fitting compared with the 
CBI2 restored MEM image (Fig. [HJ . We scaled the CBI2- 
simulated visibilities by the cross-correlation slopes obtained 
in Sec. 4.2 (see Table |3| to bring them to the same range 
as the CBI2 data, and averaged 90 reconstructions, corre- 
sponding to different realisation of Gaussian noise as given 
by the CBI2 visibility weights. This noise-simulation process 
applied to the CBI2-simulated IRAC visibilities, and subse- 
quent averaging, avoids differences in convergence with the 
models fit to CBI2, as could arise from the varying visibility 
weights. 

We see in Fig. [6] that emission from NGC 2064 is at 
24 ± 1% of that found for NGC 2068 in the CBI2-simulated 
8 /im image, while in CBI2 this percentage is 70 ± 7%. We 
also note that the radio emission peaks are not coincident 
with illuminating stars, particularly in NGC 2064 where the 
radio nebula is offset to the southeast of the infrared coun- 
terpart. If the emission is due to spinning dust, and if the 
spinning dust emissivities are fairly independent of the radi- 



ation field (e.g. |Ali-Hai'moud et al.|2009||Ysard fe Verstraete 
2009 1 , we expect the ratio of mid-IR and CBI images to be 



proportional to the averaged UV intensity. So we do not 
expect a direct correspondence between the radio and the 
IR. For instance, in the vicinity of the brightest stars the 
mid-IR intensities will be enhanced relative to the spinning 
dust intensities. The UV field will also affect the grain size 
distribution, directly impacting the spinning dust spectrum. 
Mid-IR specific intensities from dust grains are expected to 
be proportional to the UV field, as well as to very small grain 
column densities. This may explain the fact that NGC 2067 
is fainter at 8 /mi than at 31 GHz because this nebula does 
not contain stars that are bright enough to have a strong 
UV field, but still has enough dust column to radiate con- 
spicuously in the radio. 

We find that the morphological differences between 
IRAC 8 /xm and the 31 GHz continuum cannot be explained 
only in terms of the proportionality of IR emission and the 
UV field. Following the procedure described in Sec. |4.3[ we 
have estimated the UV field in M 78, parametrised as G - 
Fig. [6J: presents the simulation of CBI2 observations on the 
ratio of IRAC 8 /im to Go. It is apparent from Fig. [BJ: that 
the radio/IR correlation worsens in M 78 after division by 
Go • The peak in NGC 2068 is shifted to the south when com- 
pared to the observed from the UV field uncorrected IRAC 
image and the one observed with CBI2. We also observe 
that emission corresponding to NGC 2067 appears clearly 
in Fig. [6t, while in Fig. Rjk and b it is missing. 



4.2 CBI - IR correlations 

For a quantitative analysis we performed linear correlations 
between CBI and IRIS data and between CBI2 and IRAC 
at 8 /im, both in the visibility and sky planes. We used 
different diagnostics in each case. On the visibility plane 




Figure 6. (a) Comparison between IRAC 8 /im (in grey scale) and the CBI2 MEM restored image (in contours). The intensities are 
in MJy sr -1 (b) In gray scale we have the CBI2-simulated IRAC 8 /im image, scaled by the correlation slope with CBI2, against the 
restored image from CBI2, both reconstructed using model- fitting, (c) As in (b), but for simulations of CBI2 observations on the ratio 
of IRAC 8 fim and the UV field Go ■ Coordinates follow from Fig. [2] 



we used reduced % 2 an d the linear correlation coefficient r, 
which we calculated using: 



V N E (v 3 \) - (E vk) 2 V N E ( y t 2 ) - (E V *Y 



(i) 



where N is the number of visibilities, V31 are the visibilities 
at 31 GHz, Vt are the visibilities from the simulated tem- 
plate and the sums extend over every visibility. This last 
coefficient will be referred to as r v i s - 

The sky plane correlations were carried out within cir- 
cular areas, with radii of 45 arcmin (i.e. with diameters twice 
the CBI primary beam F WHM) , and centred on each region 
of interest (i.e. M 78 for CBI1 and NGC 2068, NGC 2071 for 
CBI2). These sky plane correlations were made with mod- 
elfitted images. For these correlations we also use the r coef- 
ficient, changing now visibilities for intensities and for the 
number of pixels. This coefficient will be referred to as r s ky 
We also use the normalized cross-correlation coefficient, 



C 



E('ii) 



(2) 



where 731 and It are the intensities at 31 GHz and from 
the template respectively and a is defined by /31 = al t and 
determined by linear regression (this definition of C is as 
used in |Casassus et al. 2006, except there is a typo in their 
Eq. 1, where a is missing). Here again, the sums are taken 
over all pixels falling within a circular aperture 90 arcmin 
in diameter. This choice of aperture is meant to include all 
of the signal seen by CBI. 

To obtain statistics on r s k y and C we averaged 1000 val- 
ues of r s k y and C taken on different reconstructions of the 
CBI data, each with the addition of different realisations 
of Gaussian noise to the visibilities (given by the visibility 
weights) . The values shown in Table [2] and Table [3] corre- 
spond to the mean value of a normal distribution fitted to 
the data so obtained, and the errors correspond to 1 a. 



Table 2. Linear Correlations Results for CBI versus IRIS, both 
for visibility and sky planes. 



Template 


X 2 /" a 


r b - 

VIS 


C c 


sky 


12 (im 


1.40 


0.776 


0.683 ± 0.024 


0.786 ± 0.019 


12 fim e 


1.39 


0.779 


0.774 ± 0.025 


0.786 ±0.018 


25 (im 


1.59 


0.768 


0.644 ± 0.023 


0.739 ±0.021 


25 fim e 


2.15 


0.719 


0.627 ± 0.022 


0.717 ± 0.017 


60 /im 


1.58 


0.755 


0.651 ± 0.020 


0.787 ± 0.017 


60 fim e 


1.28 


0.789 


0.750 ±0.024 


0.819 ±0.018 


100 ixra 


1.39 


0.777 


0.699 ±0.021 


0.802 ±0.019 


100 /im e 


1.58 


0.754 


0.761 ± 0.027 


0.761 ± 0.018 



a Reduced x 2 , taken in the visibility plane. 

v is the number of degrees of freedom, i.e. twice the number of 

visibilities (counting real and imaginary parts) 

b Linear correlation coefficient from the visibility correlations . 

c Normalized cross-correlation, from the sky plane correlation. 

d Linear correlation coefficient, obtained from the sky plane 

correlations. 

e Comparisons made with UV field corrected IRIS images. 



Table 3. Linear Correlations Results for CBI2 versus IRAC 8 (im 



Template 


x 2/ u * 


r ■ b 

' VIS 


C c 




'^sky 


NGC 2068 


3.44 


0.254 


0.473 


± 


0.017 


0.739 ± 0.020 


NGC 2068 e 


3.92 


0.129 


0.409 


± 


0.015 


0.677 ± 0.019 


NGC 2071 


3.99 


0.205 


0.347 


± 


0.028 


0.530 ± 0.034 


NGC 2071 e 


4.08 


0.180 


0.402 


± 


0.031 


0.516 ±0.036 



"Reduced \ 2 of the visibility correlation, where v is the number 
of degrees of freedom 'Linear correlation coefficient, obtained 
from the visibility correlation. c Normalized cross-correlation, 
from the sky plane correlation. ^Linear correlation coefficient, 
obtained from the sky plane correlations. e Comparisons made 
with UV field corrected IRAC images. 



Cm-wavelength continuum in M 78 9 



Table [2] summarizes numerical results from the correla- 
tions between CBI and IRIS. As we can see the best results 
independently of the diagnostic (reduced % 2 > C, r v i s or r a k y ) 
correspond to 12 /im and 100 /im, which both, given the 
uncertainties, correlate equally well with the CBI data. The 
worst correlations are obtained at 25 pm, except for the r v i s 
diagnostic, where it gives a value closer to one than the ob- 
tained result at 60 /im. 

For correlations with the finer resolution data of CBI2 
we used IRAC 8 /im as a template corresponding roughly 
to IRIS 12 /im. However the two channels do not probe ex- 
actly the same medium, we have different widths and line 
contributions. In the case of 12 fim the channel ranges from 
A = 7 /im to A = 15.5 /im (AA = 8.5 /im) and includes 
PAHs features at 7.7 and 11.3 /im and H2 lines at 8.026, 
9.665 and 12.3 /im. For 8 /im the range is from A = 6.5 /im 
to A = 9.5 /im (AA = 3 |i,m) with contributions from the 
PAH feature at 7.7 /im and H2 lines at 6.91 and 8.026 /im. 

The correlation diagnostics worsen when testing with 
CBI2-simulated data from IRAC 8 /im, as summarised in 
Table [3] The degradat ion of the correlation when comparing 
CBI2 and IRAC 8 /im might be due to the intrinsic mor- 
phological differences pointed out in Sec. |4.1| which stand 
out with the finer resolution of CBI2 and IRAC, compared 
to the CBI and IRIS tests. 



4.3 Stellar UV field and VSG column density map 

Under the spinning dust hypothesis we expect both the mid- 
IR and radio emission to be proportional to the column den- 
sity of VSGs. We assume that mid-IR emission is due to 
stochastically heated VSGs. Thus the mid-IR emission will 
also be proportional to the UV field. We also assume that 
spinning dust is proportional to the column of VSGs, and is 
fairly independent of the local UV field. 

Then to obtain a more accurate template of the VSG 
column we must divide the mid-IR images by an estimate 
of the UV field. We estimated the dust temperature Td with 
the IRIS images at 60 and 100 /im. The resulting Td map 
shown on Fig. [7]i was obtained by fitting a modified black- 
body to each pixel using an emissivity index for big grains 
of p = 2 (representative of the emissivity index in the far- 
IR) . The Td map can be converted to an equivalent radiation 
field using Go = (T d /17.5 K) /3+i (as in |Ysard et "aT||2009 1. 
The result is shown in Fig. [Tfj. 

We also fitted a modified black-body to the integrated 
flux at 100 and 60 /im obtaining an average temperature of 
32 K and a spectral index of 1.83, consistent with big grains. 
The spinning dust emission should arise from the VSG pop- 
ulation that is traced by mid-IR emission. Calculating the 
expected emission from this modified black-body at 12 /im 
we found that it is negligible (~ 10 
extracted from the 12 /im image is ' 
pect a significant population of VSGs. This is also deduced 
from the 8 /im image which shows significant emission from 
the 7.6 /tm PAH band. 

The UV luminosities of HD38563N, HD38563S, 
HD38563C, HDE290861 and SSC V 111 are extra cted from 
the ATLAS 12 model atmospheres ( |Castelli||2005[ ). The val- 
ues we used for the stellar parameters are listed on Ta- 
ble 21 In connection with the discussion on C I continuum in 
Sec. |4.5[ we include the Stromgren radius of the C n region 
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Figure 7. Proxies for the temperature and UV field in M 78. a) 
Dust temperature map obtained from IRIS 60/100, b) implied Go 
map with a dust emissivity = 2. c) IRAC 8/im filtered to remove 
point sources (i.e. stars), and resampled at IRIS resolutions, d) 
Expected Go map from the stellar content of M 78, and in the 
absence of extinction. The grey-scale has been chosen to match 
b); Go diverges as the inverse-squared projected distance to the 
UV sources. Coordinates follow from Fig. [2] 



around each star using values of n e = 1 cm and T e = 50 
K, and recombination coefficients from |Nahar fc Pradhan| 
(p97|. 



Knowledge of the stellar content can be used to estimate 
the circumstellar UV field, by diluting the stellar UV lumi- 
nosities according to the inverse-square of the projected dis- 
tance. This rough estimate of the Go map, shown in Fig. [Tjl, 
assumes that the layers of M 78 exposed to the stellar UV 
field do not overshadow. In other words, Fig. [7jl assumes 
that the mid-IR emission stems from regions directly ex- 
posed to the radiation from the illuminating stars, without 
intervening extinction. The similarity of the maps in Fig. [7)3 
and Fig. [TJi helps as a sanity check on Fig. [Tfj. 

With an estimate of the UV field in M 78, as 
parametrised by Go, we can divide the IRAC 8/im image 
by the Go map to produce a VSG column density map. 
This estimative column density map, in arbitrary units, is 
shown on Fig. [TJ;. The spike near the origin of coordinates 
in Fig. [7£ corresponds to HD38563N. 



4.4 CBI and CBI2 correlations with the VSG 
column density map 

In Fig. [8^ we see the gaussian distributions fitted to the his- 
tograms of the r s k y coefficient for the different IRIS bands 
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Table 4. Parameters for the brightest stars embedded within M 78. 



Name Spectral Type Temperature™ Surface Gravity 6 Luminosity c C II Stromgren Radius' 

HD38563N B2m 20000 -1 1.6 10 4 6.06 

HDE290861 B2v 22000 -0.5 6.6 10 3 4.84 

HD38563S B5v 15000 -0.5 6.9 10 2 1.39 

HD38563C AOn 10000 -1.5 1.3 10 3 0.11 

SSCV 111 B3v 19000 -0.5 2.7 10 3 3.16 



"In K. 

b In log(g/ S0 ). 
c InL . 

d In pc using n e = n c + = 1 cm -3 . 



before performing the UV field correction. The two gaus- 
sians fitted to the histograms at 12 and 60 /im are almost 
coincident while at 100 /im we find the best correlation, 
r s k y = 0.802 ± 0.019. Finally at 25 /im we find the worst 
correlation result, as was expected because this is the only 
template where NGC 2071 appears brighter than NGC 2068. 

In Table|2]we give the results for IRIS-CBI correlations, 
corrected by our estimate of the UV field. We see that the 
best correlation now corresponds to 60 /J,m, which improves 
from the uncorrected version, as reflected in the 3 a increase 
in r v i s , C and r s k y . For 12 /im the improvement is slight. In 
general we see that the values found for 25 and 100 /im suffer 
a worsening when compared with the UV-field uncorrected 
results. In the case of 25 /im, we see again the worst cor- 
relation results, as expected for here the brightest nebula 
is NGC 2071 and not NGC 2068. The Gaussian fit to the 
distribution of the r s k y coefficient is shown in Fig. |8j where 
we can compare cross-correlations with templates corrected 
and uncorrected by the UV field. 

The UV field corrected results for IRAC-CBI2 correla- 
tions are given in Table|3] In Fig. [9] we show the distribution 
of the correlation coefficient r s k y obtained for each nebula, 
with and without the UV field correction. In NGC 2071 
the correlation indicators are generally slightly lower after 
division by Go, but the decrement is within the noise. In 
NGC 2068 all correlation indicators decrease markedly, with 
significant decrements in C and r s ^ y over 3 a. 

The estimates of the VSG column density maps pro- 
duced by dividing the IR images by the UV-field should be 
taken with some caution. Apart from the uncertainties in 
the UV-field estimates, there is also the possibility that the 
grain population may be influenced by the UV field, in par- 
ticular in its size distribution. To make a more complete 
comparison is beyond the scope of this paper. 



4.5 Expected C I continuum 

Estimates of the C I continuum level can be made under the 
assumption of ionisation balance for the stellar content of 
M 78 with T e ranging from 20 to 50 K (see Sec. If the 
stellar U\Q radiation is absorbed by carbon in the M 78 
PDRs, and the kinetic temperature is T e = 50 K, we expect 
a flux density of F v = 2.7 Jy at 31 GHz for the whole of 



CBI/IRIS Correlation 



M 78. Details are given in Casassus et al. (2008 note there 
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Figure 8. (a) Fitted distributions to Monte-Carlo noise simula- 
tion of the correlation coefficient rgw between CBI and IRIS for 
M 78. The solid lines corresponds to the 12 fim templates, the 
dashed lines to 25 (im, the dash-dotted lines to 60 £im and the 
dotted lines to 100 fim. The Gaussian parameters are listed in 
Table [2] (b) Distributions obtained from UV-field corrected tem- 
plates. The corresponding parameters are also listed in Table [2] 



is a typo in their Eq. 1, for a missing D 2 in the denominator). 
We use the recombination coefficients tabulated in N ahar fcl 



Pradhan ( 1997 \ 



However, Casassus et al. (2008) overlooked the com- 



petition of dust and carbon for the carbon-ionising con- 
tinuum, which may be estimated by comparing opacities 
of dust, Td, and carbon, re, to C-ionising radiation. We 
take a constant photoionisation cross-section for carbon of 



( Bohlin et al. 



1978 



tion A v = N s /2 10 21 
find that r d /r c w A X /(A V 3.2) 
standard Rv — 3.1, and using the extinction curves given in 



in CGS units). We 
2.0 for A w lOOOA, for a 



the total number of carbon ionizing photons is 6.33 10 s 



Draine (20031. About one third of all carbon-ionising pho- 
tons is absorbed by carbon, although this result is sensitive 
to the exact value for Ry. 

The expected C I continuum levels at 31 GHz, in a dust- 
free nebula, are F v = 2.1, 2.7 and 3.3 Jy, for T e = 20, 50 
and 100 K. But taking into account dust absorption of the 
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Figure 9. Distribution obtained from Monte-Carlo noise simula- 
tions of the correlation coefficient r s jt y between CB12 and 1RAC 
for NGC 2068 and NGC 2071. The solid line histogram corre- 
sponds to the data obtained without the UV field correlation, 
while the dashed line one is the corrected distribution. We also 
show gaussian fits for each case, whose parameters are listed in 
Tabled 



carbon-ionising radiation, these values are reduced by 1/3, 
to not more than 1 Jy. This, however, is not negligible com- 
pared to the observed level of anomalous emission calculated 
from the differenced data set, of 2.79 ± 0.81 Jy. 

For a uniform slab nebula, seen face-on as a disk with 
a radius of 0.1 deg, and at a distance of 400 pc, the corre- 
sponding electron density at T e = 50 K is 15 cm -3 , which 
is nn — 1.5 10 5 cm -3 if the carbon in the PDR gas is all 
singly ionised. Ionisation balance indicates that the carbon 
continuum could reach observable levels, given the observed 
physical conditions in M 78 (Sec.[2|. However, the contribu- 
tion of optically thin carbon continuum at 5 GHz is probably 
negligible since the 5-31 GHz index is a 31 > (we do not 
see the nebulae in the PMN maps). To reach unit opacity at 
5 GHz, we require electron densities of order ~ 100 cm -3 , 
corresponding to H-nucleus densities of ~ 10 6 cm -3 if all of 
carbon is ionised, with a filling factor ~ 10 -4 . 

Alternatively, the lack of detectable optically thin Cl 
continuum at 5 GHz could be used to place constraints on 
the dust extinction law in M 78. This requires full-blown 
PDR models, and is beyond the scope of this work. 



4.6 Free-free specific intensities 

In Fig. [10] we compare the CBI contours and the free-free 
tracers on small scales. The 4.85 GHz PMN image traces 
Hi free-free. No 4.85 GHz emission is evident at the loca- 
tion of NGC 2071. The 4.85 GHz intensity peak, located in 
NGC 2068, is 0.073T0.011 Jy beam- 1 . If we assume an elec- 
tron temperature of T e = 7000 K, then the spectral index of 
optically thin free-free emission between 5 GHz and 31 GHz 
is a = -0.12 (e.g. [Wilson et al.|2009| pickinson et al.|2003 



Casassus et al.|2l')07;. and the peak free-free intensity should 
be ~ 0.055 Jy beam -1 at 31 GHz. Ha, also a tracer of Hi 
free-free, bears similar properties as the PMN map - no Ha 
is seen in NGC 2071, albeit faint wings from Barnard's loop. 

In each template, i.e. NVSS, SHASSA and PMN, there 
is compact free-free emission coincident with the brightest 
star embedded in NGC 2068. It is the compact Hn region 
reported by|Matsakis et aT (19761. We see that the peak is 



coincident with HD38563N for both CBI and PMN images. 

In order to make a proper comparison between PMN 
and CBI/CBI2 specific intensities we have to smooth them 
to a common resolution. Any high-pass filtering from PMN 
will not be significant in this comparison since, as explained 
in Condon et al. (19941, the PMN survey filters-out scales 
above 20 , and since the shortest CBI2 baseline is on 20 ar- 
cmin scales, CBI2 and PMN have fairly well-matched beams. 
Being an interferometer the CBI filters on all scales above its 
PSF (~8 arcmin), so we expect the comparison with PMN 
to give conservative values in terms of the detection of an 
excess over optically thin free-free. 

The specific intensity units of the PMN data retriev- 
able from Sky VievvQ are not defined precisely. It is stated 
in the header of the PMN images that 200000 beams are 
about one steradian (so a beam of 7.2 arcmin). To obtain 
the exact size of the PMN beam, and to cross check on the 
Jy beam - units, we fitted a gaussian to 3C273, which is at 
a declination similar to that of M 78 (RA 12 h 29 m 06".7 DEC 
+02°03'08".6) so both images were obtained using similar 
beams. Although an elliptical gaussian fit to 3C273 has a mi- 
nor axis of 3.8 arcmin FWHM, trials on neighbouring and 
fainter radio point sources led to values around 3.4 arcmin. 
We assume the original PMN map we are using has a reso- 
lution of 3.7 arcmin. Given the resolution of the CBI maps 
used here, of ~ 5'. 8, we smoothed the resolution of PMN to 
that of CBI by convolving with a gaussian whose FWHM 
is FWHMker = FWHM^ BI - FWHMi, MN . We also use this 
value for the beam units in the PMN image. 

After smoothing, and scaling Jy beam -1 units to the 
CBI beam, we find that the peak intensity in the M 78 
PMN image is of 0.015 ± 0.004 Jy beam -1 , while the peak 
found in the CBI image is of 0.492 ± 0.031 Jy beam -1 . Us- 
ing a spectral index a = —0.12 and assuming pure free-free 
contribution at 4.85 GHz we find an expected emission of 
0.012 ±0.003 Jy beam -1 at 31 GHz. The difference between 
the actual emission detected at 31 GHz and the one expected 
from 4.85 GHz free-free is of 0.480 ± 0.031 Jy beam -1 . This 
corresponds to a spectral index a = 1.89 ± 0.15 between 
4.85 GHz and 31 GHz. 

From archival fluxes from PMN and NVSS we find that 
the peak specific intensity between 1.4 and 4.85 GHz has a 
spectral index of 0.00 ± 0.23. This is an upper limit since 
NVSS filters out large scales (NVSS has more flux loss than 
PMN) . We can nonetheless conclude that the free-free emis- 
sion seen at 1.4 and 4.85 GHz is optically thin. 



4.7 Low Frequency Spectrum of NGC 2071 

Under the hypothesis that the absence of NGC2071 at 5 GHz 
is due to an optically thick spectrum, with the observed 
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Figure 10. CBI contours over free-free tracer images. The crosses indicate the position of the brightest stars, as in Fig. [4] We find that 
the peaks coincide, but there is no further correlation between CBI and any of the templates, (a) CBI contours against NVSS image, 
(b) SHASSA image with CBI contours. In this image we can appreciate the extended H a emission that surrounds M 78. (c) A zoomed 
image from PMN with CBI contours. Coordinates follow from Fig. [2] 



31 GHz flux, of ~0.1 Jy we expect 16 mjy beam -1 at 5GHz. 
This is just barely at 2 a in PMN, so could be missed. Also 
there is hot dust in NGC 2071, as seen in the 25 fim image, 
perhaps due to classical grains heated by a hot star. 

For optically thick emission, the expected flux density 
at 5 GHz from NGC 2071, 16 mjy beam -1 , corresponds to a 
diameter of 0.05 arcsec. Such a compact and optically thick 
source would stand out in VLA observations at frequencies 
higher than 5 GHz. Unfortunately we found no VLA obser- 
vations that cover the 31 GHz centroid of NGC 2071. 

However, the bulk of the 31 GHz flux density seen to- 
wards NGC 2071 is probably not due to an UCHII region. 
The MEM model of the CBI2 visibilities finds structure in 
NGC 2071, so that its extension is probably closer to 4'. 
If opacity is 1 at 31 GHz (so that the 5-31 GHz index is 
+2), then a distance of 400 pc and an angular extension of 
4', or a linear depth of 0.5 pc, gives an electron density of 
10 s cm -3 if T e = 10 4 K. Such densities are found only in 
the most compact UCHII regions, under 0.01 pc in diam- 
eter. In fact no UCHII larger than 0.1 pc are tabulated in 



Wood & Churchwell 



( 1989 1, in their Table 16. More recently, 



Murphy et al.| ( |2010 1 have listed UCHII region properties in 
their Table 3, where an upper limit of 0.1 pc can be found on 
the linear sizes of regions with N e w 10 s cm -3 . The ionised 
mass for NGC 2071 would have to be around 125 Mq, while 
the largest UCHII regions reach only 0.5 Mq. 



4.8 Radio Stars 

There are several protostars in this region, which contribute 
to the total emission at 31 GHz. In most cases the clumps 
associated with stellar formation do not have a related ra- 
dio continuum source (to a limit of 0.1 mjy at 8.4 GHz 



Gibb 1999). Table 2 in Gibb (19991 gives information on 



the known 8.4 GHz radio continuum sources within M 78. 
We also included a source found in the VLA archive not 



included in Gibb (19991; this sources is roughly coincident 



with NGC 2071 IRS 1. 

The radio continuum sources will contribute to some 
extent to the emission found at 31 GHz. To derive this con- 



tribution we use the formula S v oc v a , with a being the 
spectral index. Collimated, ionized stellar winds (sources for 
radio emission from young stellar objects) have spectral in- 
deces constrained by —0.1 < a < 2, with the highest val- 
ues corresponding to totally opaque sources (e.g. |ReynoIds| 
19861. Using this and considering the highest possible spec- 



tral index we find that the combined contribution from these 
radio continuum sources is 237. 98± 11.22 mjy, far below the 
observed emission at 31 GHz which is of 2.73±0.28 Jy. Thus 
we neglect any contribution from radio stars to the 31 GHz 
flux density in a 45 arcmin aperture. 



4.9 Spectral energy distribution fits 

In the previous section we used specific intensity (surface 
brightness) to make the comparison between 31 GHz and 
free-free templates. Here we use the integrated flux densi- 
ties to build a SED. For this we use the values shown in Ta- 
ble [Tjfor differenced fluxes, unless otherwise indicated (See 
Sec. |3.3| |. We will fit a model of three components of emis- 
sion: free-free, vibrational dust and spinning dust. 

In Fig |ll| we see a fit to the data with two emission 
mechanisms, not including spinning dust emission. One is 
free-free emission for which we used a = —0.12, adjusted the 
datum at 2.326 GHz. We also fitted a modified black-body 
to the IRIS data points at 100 fj,m and 60 /im. This fit gives 
us a warm dust temperature of 32 K and a spectral index 
fi = 1.83 (we take the spectral index from S v oc B V {T)). 
To check this value we also fitted a modified black-body to 
the IRIS images that were CBI-simulated, and found con- 
sistent values of T — 32 K and f3 — 1.82. We also fitted a 
Rayleigh- Jeans tail to the points of WMAP which accounts 
for cold dust, assuming that there is no contribution from 
spinning dust. The characteristic values for the spectral in- 
dex /3 (defined by S v oc v 2+fl ) of ISM dust in the Rayleigh - 
Jeans regime, ranges from 1 < P < 2 (j Shetty et al.|[20 09b). 
To measure the confidence level of our fit we used a reduced 
X 2 test. For this fit we used a value of j3 = 0.6, lower than 
characteristic values found in the ISM, but also the highest 
value for which the confidence level is ^ 5%. We see from 
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Figure 11. M 78 SED and model components. The spectrum 
of the radio flux in a 45 arcmin aperture is modeled using 2 
components made up from free-free emission (dotted line) and 
a Rayleigh- Jeans tail fitted to the WMAP data (dashed line). 
Diamonds correspond to flux densities measured by WMAP, the 
triangle is the flux density obtained by Rhode/HartRAO and the 
asterisk is the flux density from CBI-IRIS 100 (im correlated emis- 
sion. The solid line is the sum of the three emission components 
considered. 



this fit that a model which uses only Rayleigh- Jeans from 
cold dust is unable to explain the emission from M 78 at 
95% confidence. 

We us ed the spinning dust model from|Draine fc Lazar-| 



(1998) for dark clouds, which assumes nu = 10 cm 



and Td = 10 K. We fitted this model to the flux densi- 
ties from differentiated WMAP data, except the W band, 
see Table [I] For a given spinning dust emissivity per H- 
nucleus (which depends on the physical environment), the 
spinning dust intensities are determined by the column den- 
sity of hydrogen nuclei, AT H . We fit the WMAP SED by 
varying N H , and obtain a value of (7.99 ± 2.34) 10 21 cm~ 2 
for A'h = / * nu * L, where / is the filling factor and L is the 
nebular depth. The uncertainties on Nu take into account 
uncertainties on the level of free-free and the level of residual 
Rayleigh-Jeans dust emission, at 31 GHz. We see that the 
cm-wave excess, parametrised by the spinning dust model, 
is significant at (Nu) = 3Aa(Nu)- 

The value for / can be constrained using the densities 
assumed by the spinning dust model, and those given by 



Lada et al. ( 19971; Strom et al. ( 1975 I. We assume that nu 
ranges from 10 J to 10° cm _J . We then consider the neb- 
ula as a sphere, with characteristic size of 27' which is the 
angular distance from the northeast extreme NGC 2071 to 
the southwest end of NGC 2064. When we consider a dis- 
tance of 400 pc, we see that Nu would range from 10 22 to 
10 25 cm" 2 , which constrains -3.25 < log(/) < 0.01. When 
using the extinction maps from Schlcgcl, Finkbciner & Davis 
( 1998 1 we find H-nucleus column densities within a smaller 
range, with Nu = (2.28 ± 1.40) 10 23 . For this we consid- 
ered a constant ratio Ay/E(B — V) = 5.5 and th e rela- 
Av/5.3 cm" 



(Kim & Martin 19961. This 



tion Nu. = W 

new range in H-nucleus column densities gives us the value 
log(/) = -1.46 ±0.30. 

Fig.[l2]shows the spectral energy distribution calculated 
from WMAP, Rhodes/HartRAO and CBI, this last one ex- 
tracted from IRIS 100 /im using the coefficient obtained from 
the visibility plane correlation (see Sec. |3.3| l . The free-free 
emission was derived from the datum at 2.326 GHz using 



Figure 12. M 78 SED and model components. The spectum of 
the radio flux in a 45 arcmin aperture is modeled using 3 com- 
ponents made up from the vibrational black body emission of 
classical dust grains, and that from free- free and spinning dust. 
Diamonds correspond to flux densities measured by WMAP, the 
triangle is the flux density obtained by Rhodes/HartRAO and 
the asterisk is the flux density from CBI-IRIS 100 fim correlated 
emission. The dotted line is the free-free spectrum derived from 
2.326 GHz, the dashed line corresponds to the sub-mm dust tail, 
traced at 94 GHz and the dash-dotted line represents a spinning 
dust model. The solid line is the sum of the three emission com- 
ponents considered. 



a spectral index a — —0.12. We assume that the excess in 
WMAP at 94 GHz from the modified black-body fitted us- 
ing IRIS at 100 and 60 /im traces the sub-mm cool dust 
grey-body. We fixed to this datum point the Rayleigh-Jeans 
tail of vibrational dust using a characteristic /3 = 1.8 which 
lies within the range established by |Shetty et al . (2009b). 

We assume that spinning dust is traced by dust- 
correlated radio emission, using the IRIS 100 /J,m template, 
so that the 31 GHz emission is given by F^ighz = a^ioo/jm, 
where a = (0.092 ± 0.001) 10" 3 is the slope of the correla- 
tion in the visibility plane. If all visibilities were weighted 
equally, this would be equivalent to correcting the CBI flux 
density for flux loss. 

If we again use 2.326 GHz to trace free-free emission, 
but this time we differentiate to remove the diffuse back- 
ground, we find that the free-free level is 0.65 ± 0.35 Jy 
at 31 GHz. Compared with the value from dust-correlated 
emission at 31 GHz, which is of 2.73 ± 0.28 Jy, we have an 
excess of 2.08 ± 0.45 Jy over the expected level of optically 
thin Hi free-free. The emission at 31 GHz is not dominated 
by free-free, confirming the results from Sec. 4.5. 

Part of the Rayleigh-Jeans tail of sub-mm dust emission 
could reach 31 GHz. Fixing the level of the Rayleigh-Jeans 
flux density to the point at 93.5 GHz, with ft = 1.6 ± 0.3 
(which extends conservatively to shallow values), we expect 
0.40 ± 0.19 Jy at 31 GHz. We see that the sub-mm dust 
cannot account for the observed 100 /im-correlated flux den- 
sity at 31 GHz of 2.73 ± 0.28 Jy (close to the 33 GHz flux 
density from WMAP 33 GHz, of 2.79 ±0.81). When we con- 
sider the addition of vibrational dust emission and free-free 
emission, as previously derived, we find a 31 GHz excess of 
Fzi GHz = 1-68 ± 0.49 Jy, significant at 3.4 a. Comparing 
this value with the expected 2.73 ± 0.28 Jy we can say that 
at least 62±19%, of the observed 31 GHz emission cannot be 
explained by Hi free-free and/or vibrational dust emission. 

The addition of a spinning dust component improves 
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Table 5. Radio-IR correlation slopes for individual clouds in 
M 78. 



Source 



(UK (MJy sr- 1 )- 



NGC 2068 
NGC 2071 
NGC 2067 
NGC 2064 
6 H II regions" 
15 Cool dust regions 6 



4.29 ±0.25 
4.51 ± 0.14 
7.21 ± 1.72 
7.69 ±0.96 
3.3 ± 1.7 
11.2 ±5.0 



"Dickinson et al. I 2007 1. ( Davies et al. I 2006 1 



significantly the confidence level (as derived from a reduced 
X 2 test) for the SED fit, when compared with models that 
only consider the Rayleigh- Jeans tail from sub-mm dust to 
account for the spectra at 20-60 GHz. The only free parame- 
ter for a spinning dust component is f*N(H), whose value of 
(7.99 ± 2.34) 10 21 cm~ 2 differs significantly from zero. From 
this we can say that spinning dust is a significant component 
in our SED for M 78 at angular scales of ~ 8'. 



4.10 Comparison with other radio 
proportionality coefficients 



IR 



For comparison with previous work on the radio-fR corre- 
lations in other objects, we converted our correlation slopes 
to T 3 i ghz//ioo pm, i.e. fiK (MJy sr" 1 )" 1 instead of di- 
mensionless unit^3 These values are summarized in Ta- 
ble [5] We found that the highest slopes correspond to the 
fainter nebulae, NGC 2064 and NGC 2067, with values 
of 7.69 ± 0.96 and 7.21 ± 1.72 /jK (MJy sr" 1 )" 1 , while 
for NGC 2071 and NGC 2068 we found 4.51 ± 0.14 and 
4.29±0.25 fiK (MJy sr" 1 )" 1 . We see that all these values are 
within the range constrained by 3.3± 1.7 fiK (MJy sr -1 )" 1 , 
derived by |Dickinson et al.| ( |2007| from 6 Hn regions, and 
the value of 11.2 ± 5.0 (MJy sr" 1 )" 1 for cool dust ob- 



tained from 15 regions by Davies et al. (20061. Regions with 



higher temperatures tend to have lower slopes. The radio- 
correlated dust emission stems from hotter dust in Hn re- 
gions than in M 78. Also we see that the correlation slopes in 
NGC 2071 and NGC 2068 are lower than in NGC 2067 and 
NGC 2064, as expected since the exciting stars in NGC 2071 
and NGC 2068 are hotter. 



5 CONCLUSIONS 

We have found evidence of dust-correlated emission towards 
M 78 at a frequency of 31 GHz. Its spatial distribution is 
coincident with that from dust, especially as traced at 12 
and 100 fim., for which we found similar linear correlation 
coefficient in the visibility plane. The morphology at 31 GHz 
is qualitatively inconsistent with pure optically thin free- free 
emission, as can be traced by PMN, NVSS and SHASSA, 
where we only find significant emission around HD38563N. 

The radio-IR correlations worsen at finer angular reso- 
lutions. The CBI2 and IRAC 8 /im correlation tests worsen 



5 the radio-IR correlation slopes in units of (MJy sr 1 ) 1 
are sometimes referred to as 'emissivities' 



relative to those found from CBI and IRIS. The correlation 
tests further degrade when dividing the IR images by the 
UV-field. Thus the radio/IR differences cannot be explained 
as being due to modulation by the UV-field. 

The specific intensity levels seen by CBI at 31 GHz 
are higher than those expected for free-free emission de- 
rived from PMN data (smoothed to the CBI resolution) by 
0.480 ± 0.031 Jy beam" 1 . However, the flux densities mea- 
sured in a 45 arcmin circular aperture are dominated by dif- 
fuse free-free emission, as suggested by flux measurements 
with WMAP and Rhodes/HartRAO. When we differentiate 
against a reference field within the Orion B complex, we 
find an excess flux density at 31 GHz of 2.08 ± 0.52 Jy over 
the expected free-free in a 45 arcmin circular aperture. We 
need high resolution data to subtract more accurately the 
background free-free emission and confirm the values found 
in this work. 

When fitting a Rayleigh- Jeans law to WMAP differ- 
entiated fluxes we derive a spectral index /3 = 0.6. This 
value is lower than those typically found in the ISM, where 
1 < P < 2, but even allowing for this value of /3 the confi- 
dence level of this model spectrum, given the observed SED, 
is only of 5 %. We can neglect the contribution of a Rayleigh- 
Jeans tail from sub-mm dust at cm-wavelengths. 

The addition of a spining dust component to the SED 
fit results in a significant improvement in confidence level. 
The only free parameter in our fit, the column of H- nuclei 
which scales the model spinning dust emissivities per H- 
nucleus, differs significantly from zero, thus confirming that 
the higher confidence level achieved is due to a non neg- 
ligible component. This is also supported by the fact that 
the combined contributions of Rayleigh- Jeans and free-free 
at 31 GHz yield a total flux that is 1.68 ± 0.49 Jy below 
the emission detected by CBI. A spinning-dust-like compo- 
nent is needed in order to explain the observed emission at 
cm- wavelengths in M 78. 
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